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Abstract 

We adopt a Markov Chain Monte Carlo method to examine various new physics models which can 
generate the forward-backward asymmetry in top quark pair production observed at the Tevatron 



1^ , by the CDF Collaboration. We study the following new physics models: (1) exotic gluon G' , 



(2) extra Z' boson with flavor-conserving interaction, (3) extra Z' with flavor-violating u-t-Z' 
interaction, (4) extra W with flavor-violating d-t-W' interaction, and (5) extra scalars S and 



. with flavor-violating u-t-S and d-t-S^ interactions. After combining the forward-backward 

I asymmetry with the measurement of the top pair production cross section and the tt invariant 

cn 

' mass distribution at the Tevatron, we find that an axial vector exotic gluon G of mass about 



1 TeV or 2 TeV or a W' of mass about 2 TeV offer an improvement over the Standard Model. The 
other models considered do not fit the data significantly better than the Standard Model. We also 



?H ' emphasize a few points which have been long ignored in the literature for new physics searches: 

(1) heavy resonance width effects, (2) renormalization scale dependence, and (3) NLO corrections 
to the tt invariant mass spectrum. We argue that these three effects are crucial to test or exclude 
new physics effects in the top quark pair asymmetry. 
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I. INTRODUCTION 



The CDF Collaboration has observed a 2.3(T deviation in the forward-backward (F-B) 
asymmetry of top quark pair production at the Tevatron, using a data sample with 3.2 fb^^ 
integrated luminosity 

A%{cos 6) = 0.193 ± 0.065(stat) ± 0.024(syst). (1) 

This measurement improves the previous CDF result based on 1.9 fb~^ 

A%{cose) = 0.17 ± 0.08, and Af^i^ri) = 0.24 ± 0.14, 

where the results given in the lab (pp) and the center-of-mass (cm.) frame of the top quark 
pair (ti) are consistent with the theoretically expected dilution of 30% in passing from ti to 
pp It is also consistent with the DO result based on 0.9 fb~^ j4|: 

= 0.19 ± 0.09(stat) ± 0.02(syst), and A^^^ = 0.12 ± 0.08(stat) ± O.Ol(syst) 



for exclusive 4-jet events and inclusive 4-jet events, respectively. Although t 
consistent at a confidence level of ~ 1.5% with the SM prediction, which is 



le value is still 
6| 



A%{cose) = 0.051 ± 0.015, (2) 

it is interesting to ask whether or not the large central value can be explained by new physics 
(NP) after one takes into account other Tevatron experimental measurements of top quark 
pair production. There has been recent excitement among theorists for this measurement 



at the Tevatron |7H17|. 



In this work we point out that a strong correlation exists between ApB and a(tt) mea- 
surements and further derive the bounds on NP from both measurements under the inter- 
pretation of a variety of models. 

One should also keep in mind that, thanks to pp collisions, the Tevatron offers the best 
opportunity for measuring the asymmetry of top quark pair production, because of the basic 
asymmetry of the production process. At the Large Hadron Collider (LHC), the asymmetry 
of top quark pair production is an odd function of the pseudorapidity of the ti pair, due to the 
lack of definition of the forward direction. Hence the LHC will improve the measurement 
of the total cross section of top quark pairs, but has very limited reach for studying the 



asymmetry. In this sense the Tevatron plays a unique role for testing top quark interactions, 
and it would provide more accurate measurements with future accumulated data. Projected 
bounds on both ApB and a{tt) at the Tevatron with 10fb~^ integrated luminosity are also 
presented. 

The paper is organized as follows. In Sec. II we examine the correlation between Aps 
and a(ti) based on the recent Tevatron measurement, using the Markov Chain Monte Carlo 
method. We then give examples of a few interesting NP models generating the asymmetry, 
e.g., an exotic gluon G' (Sec. IIIip . a model-independent effective field theory approach (Sec. 
HVT) . a flavor-conserving Z' boson (Sec. |V]), a flavor- violating Z' or W (Sec. rvT|l . and a new 
scalar S{S^) (Sec. EHD. We then conclude in Sec. IVnTl 



II. CORRELATION OF ^^^b AND a{tt) 



The asymmetry ApB the top quark pair production can be parameterized as follows: 

SM _ SM , NP _ NP 
Atot ^ '^F "B + "F /q\ 

^FB SM _|_ ^SM I ^NP _|_ ^NP 

Op tOj3 -\- O p -T O £j 

NP _ „NP / _ ^SM \ ^NP 

X ( 1 4- ^ B \ (^tot 

^NP ^NP ^SM , ^NP \^) 



-^B 1 I (hot 

c^r +(^r [ - ^^r ) ^fof + <r 

\NP D , aSM 



A%xR + A'p'j^{l-R) (5) 



where 



NP _ NP SM _ ^SM ^NP 

aNP ^ ASM ^ (^F Ob „^ 1 p _ "tof /^iN 

— ^NP I ^TVP ' — ^SM I ^5Af '^^^^ ~ ^SA/ I ^TVP > 

(jp -r (Jb Op -\- atot "T o^tot 



is the asymmetry induced by the NP, the asymmetry in the SM, and the fraction of the NP 
contribution to the total cross section, respectively. In this work we consider the case that 
the NP contribution io ApB occurs in the process qq — > ti, for which the SM contributions 
do not generate any asymmetry at all at LO. However, at NLO a nonzero Af.^ is generated. 

It is worth while emphasizing the factorization of A^^ and R in Eq. ([5]), as it clearly 
reveals the effects of NP on both the asymmetry and the top quark pair production cross 
section. For example, when NP effects generate a negative forward-backward asymmetry, 
they still produce a positive observed asymmetry as long as they give rise to a negative 
contribution to cr(tt). This is important when the effects of interference between the SM QCD 
channel and the NP channel dominate. Moreover, the possibility of negative contributions 
to ap^ or a'^^ means that \Ap^\ can exceed 1. 



3 



Recently, the CDF collaboration [18| has published new results on the ti cross section 
in the lepton plus jet channels using a neural network analysis, based on an integrated 
luminosity of 4.6 fb~^, 

a{mt = 171.0 GeV) = [8.33 ± 0.40(stat) ± 0.39(sys) ± 0.17(theo)] pb , 
a{mt = 172.5 GeV) = [7.63 ± 0.37(stat) ± 0.35(sys) ± 0.15(theo)] pb , 
a{mt = 175.0 GeV) = [7.29 ± 0.35(stat) ± 0.34(sys) ± 0.14(theo)] pb , (7) 

and also an analysis combining leptonic and hadronic channels with an integrated luminosity 



of up to 4.6 fb"^ [191 



a{mt = 172.5 GeV) = [7.50 ± 0.31(stat) ± 0.34(sys) ± 0.15(theo)] pb . (8) 

Note that the theory uncertainty is derived from the ratio with respect to the Z cross section 
and the central value is quoted after reweighting to the central values of the CTEQ6.6M 
PDF 20|. By means of the ratio with respect to the Z cross section, the luminosity- 
dependence of the theoretical tt cross section is replaced with the uncertainty in the theo- 
retical Z boson production cross section. That reduces the total uncertainty to 7%, greatly 
surpassing the Tevatron Run II goal of 10%. 

In this work we fix the top quark mass to be 175 GeV as we also include the the CDF 
measurement of the invariant mass spectrum of top quark pairs in our study, which is based 
on nit = 175 GeV. We rescale the combined CDF measurements at rrit = 172.5 GeV (cf. 
Eq. [H]) to = 175 GeV which we estimate to be 

a{ti) = [7.0 ± 0.5] pb, (9) 



on the basis of the approximate 
by Langenfeld, Moch, and Uwer 



behavior of Eqs. [7] and [8] and the theoretical calculation 
2l|. It yields \R\ < 7% at the la level. Any asymmetry 
induced by the NP (Ap^) is highly suppressed by the SM cross section due to the small 
fraction R; see Eq. ([5]). 



A. Parameter estimation 



In this work we utilize a Markov Chain Monte Carlo (MCMC) to examine the correlation 
of ApB and R. The MCMC approach is based on Bayesian methods to scan over specified 
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input parameters given constraints on an output set. In Bayes' rule, the posterior probability 
of the model parameters, 6*, given the data, d, and model, M , is given by 

7r{e\M)pid\e,M) 

pmM) = — — , (10) 

where 7r{6\M) is known as the prior on the model parameters which contains information 
on the parameters before unveiling the data. The p{d\6,M) term is the likelihood and is 
given below in Eq. [TTJ The p{d\M) term is called the evidence, but is often ignored as 
the probabilities are properly normalized to sum to unity. In using the MCMC, we follow 
the Metropolis-Hastings algorithm, in which a random point, 6i, is chosen in a model's 
parameter space and has an associated likelihood, £j, based on the applied constraints. A 
collection of these points, {6i}, constructs the chain. The probability of choosing another 
point that is different than the current one is given by the ratio of their respective likelihoods: 
min(%i^, 1). Therefore, the next proposed point is chosen if the likelihood of the next point 
is higher than the current. Otherwise, the current point is repeated in the chain. The 
advantage of a MCMC approach is that in the limit of large chain length the distribution 
of points, 6i, approaches the posterior distribution of the modeling parameters given the 
constraining data. In addition, the set formed by a function of the points in the chain, 
f{6i), also follows the posterior distribution of that function of the parameters given the data. 
How well the chain matches the posterior distribution may be determined via convergence 
criteria. We follow the method outlined in Ref. 
25000 unique points in the chain. 
We adopt the likelihood 

where yij are the observables calculated from the input parameters of the z*'* chain, dj are the 
values of the experimental and theoretical constraints and aj are the associated uncertainties. 
In our case, the input parameter set is taken to be 9i = {crf^^'^ , cr^^, Ap^, Ap^}. We scan 
with flat priors for the unknown inputs with a range of 

-5 pb < af/ < 5 pb 
-5 < A^E < 5 

(recall that as a result of its definition, |A;pbI exceed 1), while the known inputs are 



22| to verify convergence after generating 
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scanned with normal distributions about their calculated central values 

a^M ^ g 3g ^ 0.50 pb 



AfJ^ = 0.051 ± 0.015 



(13) 

The calculated total ti production cross section at NLO for mt = 175 GeV has been taken 



as 
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25\ 

atH{tt) = 6.38^[!:?(scale)+S:J (PDF) pb, (14) 

where the PDF uncertainty is evaluated using the CTEQ6.6M PDF y . The fully NNLO 
QCD correction to top pair production is highly desirable to make a more reliable prediction 
on the asymmetry. Since it is still not clear how the asymmetry will be affected by the 
complete NNLO QCD corrections, we consider the NLO QCD corrections to the top quark 
pair production throughout this work without including the partial NNLO QCD corrections 



computed in 
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27|. 
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While both CDF and DO have measurements of the invariant mass distribution 
only CDF presents an unfolded differential cross section. Therefore, we inspected the ti 
invariant mass spectra reported by CDF; see Fig. [1] We take the 7 bins with Mti > 400 
GeV in our fit and weight their by the number of included bins. This assigns an equal 
weight between the Mf^ measurement and the and Aps measurements. 

The observables are di = {a(tt)'^^P, A^p^ } in addition to the binned data and define 
the output set. We use the combined cross section of Eq. ([9]). We therefore assign di = 
{7.00, 0.193} and a-i = {0.50, 0.069} in our implementation of the likelihood defined above 
for the case we denote as "Current" {J Cdt = 4.6 fb"^ for aa, J Cdt = 3.2 fb"^ for ApB, 
and J Cdt = 2.7 fb"^ for ^) while di = {7.00, 0.193} and = {0.34, 0.039} for the case 
we denote as "Projected," where / Cdt = 10 fb^^ of integrated luminosity is used for each 
measurement, in which we assume the central values are fixed and the uncertainties are 
scaled by a factor l/\/C. We combine the chains to form iso-contours of la, 2a, and 3a 
significance via their respective p-values. 

For illustration, we plot these contours in the plane of and a{tt) in Figs. EJ^a). We 
note that the current average values oi ApB and a{tt) are consistent with the SM within the 
2a level. With an upgraded integrated luminosity of / Cdt = 10 fb"^ at the Tevatron, the 
statistical uncertainty would be reduced significantly; see Fig. |2]J^b). The deviation of Aps 
from zero is then larger than 3a. Note that we also allow negative values of in this 
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Unfolded [GeV/c ] 



FIG. 1: The invariant mass spectrum measured by CDF assuming = 175 GeV. The sohd 
histogram is the CDF expectation taken from a LO calculation and Pythia. The solid band 
indicates the full NLO SM prediction with a theoretical error due to scale uncertainty which we 
use in our scans. The dashed line is K {dui^o / dMn) with K = (Tnlo/clo which shows a large 



deviation from the data. The data are taken from Ref. 



work, though they are not preferred. Taking the SM theory prediction, we translate (T(tt) 
into R defined in Eq. ([6]). The correlation of and R is shown in Figs. Ws) and (d). 
Finally, using Eq. (|5]), we obtain the correlation between and R shown in Figs. [2]^e) 
and (f). Clearly, the smaller R the larger A^^] see the la contour (solid black). 

Note that since the MCMC is sensitive to the relative likelihood change in going between 
two points, it is sensitive to only the Ax^ values. Therefore, the iso-contours of the p-values 
for 1, 2, and 3cr assume the given model. To obtain an overall indication of how well the 
model in question fits the data, we quote {x^) chain, the psr degree of freedom values 
averaged over the entire chain. In cases where we include the -^^-zlbm constraint, N^of = 3, 
otherwise Ndof = 2. This quantity is an overall estimate of the model's consistency with the 
data. Generally, values of (x^) chain < 2 are considered fairly good fits, while values much 
beyond that are not considered very good. 

One might be tempted to search for the parameter set that yields the best fit to the given 
data. However, this is doing so without regard to the level of fine-tuning required to find 
such a point. Explicitly, this can be seen as a set of points in parameter space by which 
the value is minimized, ideally to zero. However, if a small deviation from these points 
provides a large increase in x^, this particular set of points that provide a good fit can be 
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o-™{pb) R 



FIG. 2: (a) Correlation between observed A^p*^ and a(tt) at the Tevatron with the current integrated 
luminosities; (b) same for / Cdt = 10 fb^^. (c) Correlation between observed and fraction of 
NP contribution to a{tt) for the current integrated luminosities; (d) same for / Cdt = 10fb~^. (e) 
Correlation between the NP-generated asymmetry Apj^ and R for the current integrated luminosi- 
ties; (f) same for j Cdt = 10fb~^. Here, we do not include the M^f spectrum constraint. In each 
case the predicted correlations for J Cdt = lOfb"^ assume the same central values of experimental 
data. The crosses correspond to the Standard Model predictions of the Aps and cr^f. Innermost 
contours correspond to la deviation from minimum-x^ solutions; next-to-innermost correspond to 
2a; outermost correspond to 3a. Note that these scans are performed by comparing only to the 
measurements of the total cross section and asymmetry and not the M^i distribution. These should 
be compared to a SM scan only subject to these two measurements which gives {x'^)chain ~ ^'^^ 
for the current luminosity and (x^) chain ~ fb~^ if the central values are unchanged. 
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seen as more fine-tuned compared with another solution set without such a steep increase in 
X^. Therefore, the MCMC approach does take into account the parameter space available 
that affords a good fit, preferentially solutions with low fine-tuning. 

To compare the MCMC results of Fig. [2] and subsequent Figures, we ran a MCMC with 
a pure SM explanation by explicitly setting Ap^ and a^^ to zero and scanning over Eq. [T3] 
with gaussian priors. We find that 

„ f 2.56 Current Luminosity , , 

{xrjL = { (15) 

I 6.20 Projected : J Cdt = 10 fb~^ 

where "Current Luminosity" refers to the measurement of au with an integrated luminosity 
of 4.6 fb~^ and the measurement of Aps with an integrated luminosity of 3.2 fb~^. For the 
projected integrated luminosity of 10 fb~^, we assume the central values of and a^f 
remain unchanged from the values taken in Eqs. |2] and [9] while the uncertainties reduce by 
a scale factor a/£. When we examine specific models that could give rise to a larger ApB 
than the SM, we must also take into account the distribution measurement. To compare 
these models against the SM, we again run a MCMC with a pure SM explanation scanning 
over ApB and au as above while also scanning over our NLO prediction (seen in Fig. [1]) 
with gaussian priors for the last seven bins of the CDF distribution. The bin nearest tt 
threshold accounts for the majority of the total cross section. Since we already include the 
total cross section in our fit, we do not include this bin in our fit of the Mu distribution so 
that we do not weight the total cross section too heavily. If we include the measurement of 
the Mtt distribution and perform a MCMC scan over the SM, we find 

2 f 1.75 Current Luminosity 

(x )chain = \ 1 • (16) 

U.22 Projected : /m = 10 fb"^ 

Here, "Current Luminosity" refers to the above values of integrated luminosity for the au 
and ApB measurements and 2.7 fb^^ for the measurement of the Mu distribution. For the 
projected luminosity of 10 fb~^, we again assume that the central values of all measurements 
remain the same while their errors scale as l/\/C. We note that the values of {x^) chain in 
Eq. f[T^ are less than those in Eq. f[T5]l . This is because (x^) chain is a per degree of freedom. 
There are two degrees of freedom in Eq. (1151) and three in Eq. (1161) with the addition of the 
Mu distribution. The good agreement of the Mu distribution in the SM with data (seen in 
Fig. [1]) causes the per degree of freedom to decrease when it is included in the fit. When 



comparing models, we can say that if the {x^) chain value for a given model is less than that 
for the SM with the appropriate data into account, the model will provide a better overall 
fit to the data than the SM. 

The forward-backward asymmetry, defined in terms of a ratio of cross sections, is very 
sensitive to the renormalization and factorization scales, and respectively, at which 
the cross sections are evaluated. The uncertainties in the cross section associated with those 
scales can be considered as an estimate of the size of unknown higher order contributions. 
In this study, we set = = f^o and vary it around the central value of fiQ = nit, where 
rrit is the mass of the top quark. Typically, a factor of 2 is used as a rule of thumb. Large 
scale dependence in the LO cross section can be significantly improved by including the 
higher order QCD and EW corrections. In this work we calculate the SM top pair produc- 
tion cross section with the NLO QCD corrections. Unfortunately, the QCD corrections to 
the G' induced top pair production are not available yet. Therefore, we calculate the NP 
contributions only at LO and rescale them by the Mff-dependent SM qq K-factors. Due to 
the mismatch between the SM and NP cross sections, ApB calculated in this way depends 
on the choice of scale. 

Table [T] shows the LO and NLO top quark pair production cross sections in the SM at 
the Tevatron. We present the quark annihilation and gluon fusion processes individually as 



TABLE I: Predicted LO and NLO top pair production cross sections (pb) in the SM at the Tevatron, 
with = mt = 175 GeV. Note that the negative a{gq) and a^gq) cross sections are due to the 



negative NLO PDFs 



201]. 





LO 


NLO 


/tio/2 /Uo 2/io 


/io/2 /io 2/io 




6.82 5.01 3.79 


5.70 5.56 5.04 




0.37 0.24 0.17 


1.00 0.90 0.74 


(^{gq) 


0.00 0.00 0.00 


0.01 -0.03 -0.05 


o-igq) 


0.00 0.00 0.00 


0.01 -0.03 -0.05 




7.19 5.26 3.96 


6.72 6.39 5.69 


<y{gg)/(^tot 


0.05 0.05 0.04 


0.15 0.14 0.13 






0.93 1.22 1.42 
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TABLE II: Predicted LO and NLO tt production cross sections (pb) for various top quark masses 
and three scales (/io/2, fJ^o, 2/io with fiQ = mt) at the Tevatron. 





LO 


NLO 




;/n/2 


1 In 


2lln 




tin 


2;;n 


171.0 


O. WO 


5 91 


4 45 


7 61 


7 93 


6 44 


172.0 








7 *^7 


7 01 




172.5 


7.72 


5.66 


4.26 


7.25 


6.90 


6.14 


173.0 


7.61 


5.57 


4.19 


7.14 


6.79 


6.05 


174.0 


7.40 


5.42 


4.07 


6.92 


6.58 


5.86 


175.0 


7.19 


5.26 


3.96 


6.72 


6.39 


5.69 


176.0 


6.98 


5.11 


3.84 


6.51 


6.19 


5.51 


177.0 


6.78 


4.96 


3.73 


6.31 


6.01 


5.35 



well as their sum. The CTEQ6.6M 20[ and CTEQ6L 30| PDF packages are used in the 
NLO and LO calculations, respectively. In the last row we also list the K-factor, defined as 
the ratio of NLO and LO cross sections, for three scales. 

We argue that the higher order corrections cannot be estimated by a K-factor (defined as 
the ratio of NLO and LO cross sections) because the K-factor is very sensitive to the scale. 
Furthermore, the gluon fusion channel contributes much more at the NLO (roughly about 
13 ~ 15% of total cross section) than at the LO (only about 5%). Hence, one also needs to 
take account of the gluon fusion channel contribution when calculating Afb- 

Another uncertainty originates from the top quark mass. In Table Ull we show the top 
pair production cross section for various top quark masses and three scales. The central 
values of NLO theory calculations for the three masses = 171.0, 172.5, 175.0 GeV are 
always O (la) below the recent CDF results given in Eqs. ([7]) and (IHj), suggesting that the 
NP should contribute positively to ti production. 

In the following sections, we study a few interesting new physics models which can gener- 
ate a significant deviation from the SM expectation for A^r^ in the top quark pair production 
channel. We also comment on the scale dependence in each new physics model. Without 
losing generality, in the rest of this paper, we set nit = 175 GeV. 
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III. EXOTIC GLUON 



We begin with an exotic gluon (G') model, as the other models can be easily derived 
from the G' model result. In Sec. IIII Al we present analytic formulae for au and ApB- We 
calculate its width in Sec. UTTBl In Sees. UTTQ UTTDl and UTTEl we perform MCMC scans 
over parameters in several G' scenarios subject to the experimental constraints. 



A. Differential cross section and asymmetry 

The G' boson couples to the SM quarks also via the QCD strong interaction, 

G'qq : tgst'^Y {hPL + IrPr) , (17) 
G'tt : zg.t^Y {QlPl + 9rPr) , (18) 

where we normalize the interaction to the QCD coupling, Qs-, and use q to denote light quarks 
of the first two generations. Such an exotic gluon can originate from an extra-dimensional 
model such as the Randall- Sundrum (RS) model Q], chiral color model 31 40|, or top 



composite model [4]|. As discussed below, the axial coupling of G' to the SM quarks is 
necessary to create a forward-backward asymmetry. In the extra-dimensional model, such 
non-vector coupling of Kaluza-Klein gluons to fermions arises from localizing the left- and 
right-handed fermions at different locations in the extra dimension. 

The differential cross section with respect to the cosine of the top quark polar angle 9 in 



the tt center-of-mass (cm.) frame is 

= -Asm + A^i^j- + A%ps (19) 



MG') , , ,G' , AG' 



where 



dcos9 



Asm = ^{2-(3' + i(3cosef), (20) 
AG' _ ^/^«' Hs-ml,) 

2-/5') + ^ + ^)'| ' (21) 

^ [JL + JR) [9l + 9r) ) 

^ + 4^ (l - P') + 2 |^;|||;g /3cos^ + (/3cos^)^ j . (22) 



9l + 9r 



12 



Here the angle 9 is defined as the angle between the direction of motion of the top quark 
and the direction of motion of the incoming quark (e.g., the u-quark) in the ti cm. system. 
The subscripts "SM" , "INT" and "NPS" denote the contribution from the SM, the interfer- 
ence between the SM and NP, and the NP amplitude squared. For the G' model, the SM 
contribution is from the gluon-mediated s-channel diagram, the NPS contribution from the 
exotic gluon G'-mediated diagram, and the INT contribution from the interference between 



the two. The squared cm. energy of the ti system is s = {pq + Pqf , and (3 = y^l — 4mf / s 
is the top quark velocity in the ti cm. system. 

The forward-backward asymmetry of the top quark in the ti cm. frame is defined as 

A% = (23) 

where 

ap = "3 -acosti, as = -a cos 6. (24j 



/o dcosO J-i dcos9 

We further parameterize the differential cross section da /d cos 6 as follows: 

= Ai + Bicose + Qcos^e, (25) 

d cos 6 

where the subindex i denotes "SM", "INT" and "NPS". Hence, after integrating over the 
angle 6*, we obtain the asymmetry and total cross section 

^^^-^Mwy .«^i:(2^, + 2^). (26) 

where the sums are over the SM, INT and NPS terms. In reality the incoming quark could 
originate from either a proton or an anti-proton, but it predominantly comes from a proton 
due to large valence quark parton distribution functions. Taking the quark from the anti- 
proton and the anti-quark from the proton contributes less than 1% of the total ti cross 
section. Therefore, in pp collisions at the Tevatron one can choose the direction of the 
proton to define the forward direction. 

Now let us comment on a few interesting features of the asymmetry and cross section 
generated by the INT and NPS effects individually, because both effects are sensitive to 
different new physics scales: the former to a higher NP scale and the latter to a lower scale. 
First, we note that the asymmetry is sensitive to the ratio of coupling (squared) differences 
and sums for the INT (NPS) effects, e.g., 

AiNT ^ Ul- fR){9L- 9r) 2{/3) 

ih + fR)i9L + 9r) 2(2 - + 1 m ' 
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,NPs ^ ifl-PnM-gl) ^ 

ifl + pR)i9l + 9l) 2[l + il-{(3^))i2g,gn)/igl + gj,)] + l{f3r ^ ^ 

where and {(3^) denote the averaged /3 and after integration over the angle 6 and 
convolution of the partonic cross section with parton distribution functions. 

To make the physics source of the asymmetry more transparent, we define the reduced 
asymmetry {Aps) and reduced cross section a as follows: 

2INT ^ AiNT I (/l - ]r){9l - 9r) /2q\ 
'''' I ifL + fR)i9L + 9R) ^ ' 

A NFS _ A NFS 

~ / Ul + PR){91 + 91) 



a 



'INT ^ 



INT 



{fL + fR){9L + 9r) 



(31) 



NFS 

p^NFS ^ '± (on\ 

{PL + PR){91 + 91) 

The reduced asymmetries and cross sections are easily computed universal functions that 
allow us to focus on two separate limiting cases; the new physics contribution to ati and 
ApB '^^ primarily from the INT term if it is produced by a heavy resonance that interferes 
with the SM production process. If the new physics is due to a resonance that doesn't 
interfere with the SM production, then the new contribution to (Tu and Aps is given by 
the NPS term. One simply has to multiply the reduced asymmetry or cross section by 
the appropriate combination of couplings to obtain the full new physics contribution to ati 
and ApB- In Fig- [3]we plot the reduced asymmetry (a) and the reduced cross section (b) 
as functions of mc' for various choice of TG'/fnG'- The reduced asymmetry generated by 
the INT effects increases rapidly with increasing mc and finally reaches its maximal value 
~ 0.4. The reduced asymmetry generated by the NPS effects is large, typically around 
0.6-0.7. As expected, the reduced cross section of the NPS effects is alway positive; cf. the 
upper three curves in Fig. [3](b). On the other hand, the reduced cross section of the INT 
effects is always negative due to {s — rriQ,) in the numerator of Eq. fl2T|) . Both reduced cross 
sections, especially the NPS effects, are sensitive to the G' decay width. They both go to 
zero when G' decouples. 

The difference between the two reduced asymmetries can be easily understood from the 
P distribution shown in Fig. |H Figure Hl^a) shows the normalized differential cross section 
with respect to /3 for 175 GeV top quark production in the SM at the Tevatron, which 
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Mc (GeV) Mc (GeV) 



FIG. 3: (a) Reduced asymmetries defined in Eqs. (I29p and (j30p : A^pb" (lower three curves) and 
Apj^^ (upper three curves); (b) Reduced cross sections defined in Eqs. ([?T|) and ct^^-^ (lower 
three curves) and a^^'^ (upper three curves). 




FIG. 4: (a) Normalized differential cross section for production as a function of /? = y 1 ~ Amf/s 
of a 175 GeV top quark at the Tevatron; (b) (3cS = y^l — 4m^ / m^, as a function of heavy resonance 
G' mass. 
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peaks around /3 ~ 0.65. The INT effects only slightly shift the peak position. Substituting 
~ 0.65 into Eqs. ([27]) and ([29]), we obtain A^/^^^ ~ 0.4. On the contrary, the NPS effects 
prefer a much larger /3 enforced by the heavy G' resonance. We plot f3cs = \/l — 4mj / m^, 
in Fig. m^b), where ~ 0.98 for a 2 TeV G". Such a large leads to the large value 
of Apj^^ in Fig. |3l For an extremely heavy G', /3efr is equal to 1, yielding the well-known 
maximal Ap^^ = 3/4. When both INT and NPS effects contribute, one cannot factorize out 
the couplings as in Eqs. (1 27113 2 p due to the presence of both linear and quadratic coupling 
terms. 



B. G' decay width 



The An PS term contributes significantly in the vicinity of mc where the decay width 
Tg' plays an important role. Hence, it is very important to use an accurate decay width in 
the parameter scan. We consider the case that the G' boson decays entirely into SM quark 
pairs, yielding the following partial decay width 



a. 



r{G' ^ tt) = 



m 



QgLgR- 



G' 



m 



G' 



a. 



r(G'^bb) = ^[gl + gl) 



(33) 

(34) 
(35) 



where Nf = 4 denotes the light quark flavors and we have assumed that bji couples to G' 
with the same strength as t^j. In the limit of Mq' ^ mt, the total decay width of G' is 



rc' = ^ma> [{gl + gl)+2{fl + l 



(36) 



When the couplings gL/B. — Il/r are of order 1, Tg' — OigMG' — 0.1 Mg". When gi/R ~ 
fi/R ~ 3, Tg' ~ Mg'- In the following parameter scan we vary the couplings of the G' 
boson in the range of —3 to 3 when all couplings are present but —5 to 5 when only two are 
non-zero. 



C. Left-handed G': fR = gR = {) 



Since there are five independent parameters (four couplings and the G' boson mass) in 
Eqs. (I19tl22p . we turn off the right-handed couplings and gR in order to make the physics 
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FIG. 5: Theoretical prediction contours in the plane of /l and gi with fR = gR = for a 1 TeV G': 
(a) contours of cross section (pb) induced by the INT effects, (b) contours of cross section induced 
by the NPS effects, (c) contours of net cross section of the INT and NPS effects, (d) contours of 
the observed asymmetry, (e) contours of Tq/ / niQ/ . 



origin of the asymmetry more transparent. We first consider rriG' = 1000 GeV in Sec. IIII C II 
and then rriG' = 2000 GeV in Sec. IIII C 21 We will comment on non-zero /r and in 
Sees. [mD] and [mB 



1. mG' = 1000 GeV 

We first examine theoretical predictions of A°p^ and ajNT and aNPs before we perform 
a MCMC scan over the parameters. By "theoretical" we mean that the asymmetry and 
the top pair production cross section are calculated independently, without regard to their 
correlation. Fig.O^a) displays the cross section contours generated by the INT effects (ct/at't). 
The INT effects could be either positive or negative, depending on the sign of the coupling 
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product /lS'l- The INT effects dominate in the region of \fl < mG', so their contribution 
to the top pair production cross section can be written as 

a I NT oc - (m^, - sj fLQL- (37) 

This expression thus yields a positive contribution to the cross section when figh < (i.e., 
the second and fourth quadrants in Fig. [5]) and a negative contribution when figL > (i.e., 
the first and third quadrants). 

On the contrary, the NFS contribution is always positive; see Fig. [5]^b). Since the NFS 
effects contribute mainly in the vicinity of m^/, i.e. s ^ rriQ,, their contribution to the top 
pair production cross section can be written as follows: 

aNPS oc ^ ~ ml,i2fl + gir ^'^^ 

Hence, the contour pattern of cnps is determined by f1g\/ {2fl + g\Y. Fig. [5t^c) shows the 
competition between the INT and NFS effects. For a 1 TeV G', the INT generally dominates 
over the NFS. Note that the contours of net cross section are not symmetric between 
and gi due to the width effects. In Fig. |5]|^d), we see that, except for small couplings in 
the upper-right and lower-left quadrants, a positive asymmetry can be generated in all four 
quadrants. In the upper-right and lower-left quadrants, for \fL \ ,\gL\ ~ 3, the NFS term 
is large enough to generate a positive Aps- Note that negative values of ApB, although 
not plotted in Fig. [5t^d), are still consistent with the Tevatron data within 3cr C.L. and are 
included in the following analysis. 

Now we perform a MCMC scan over the parameter space after combining measurements 
of the tt asymmetry, the total cross section, and To obtain the Mu distribution, we 

separate the contribution from the qq initial state (which includes the NF that we analyze) 
from that of the gg initial state, noting that the gq and gq contributions are negligible as 
seen in Table [11 We multiply these leading order results by the SM K-factors, Ki^g and Kgg 



42l | to calculate 



respectively, which are obtained by using the Monte Carlo program MCFM 
the full NLO SM differential cross section. Each K-factor itself differs as a function of Mh fas 
seen in 



43| ) and so we weight each bin in the Mn distribution by the appropriate K-factors. 
We vary the scale /io at which we evaluate the NLO differntial cross section between mt/2 and 
2mt which gives a range of K-factors for each Mti bin. This is used in our fits as our estimate 
of the theoretical uncertainty. This uncertainty is about 10% in the first six bins and around 
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FIG. 6: Correlation of couplings for mc' = 1000 GeV with fR = gR = 0. Figures on left are 
for the current integrated luminosity while those on rig ht are for J Cdt = 10 fb"^ The regions 
inside the contours are preferred, corresponding to la (innermost), 2a (next-to- innermost), and 
3a (outermost). The crosses correspond to fL = gL = 0. 

15% in the last bin. Observe that this procedure, when NP effects are decoupled, reproduces 
the exact NLO SM differential cross section seen in Fig. [H In this and subsequent Figures of 
MCMC distributions, we adopt fiat priors in all variables scanned. The priors for the SM- 
only contribution to the ti cross section and ApB are given in Eq. f|T3l) . Fig. |6t^a) displays the 
parameter space consistent with both measurements at the la (innermost region), 2a (next- 
to-innermost region) and 3a (outermost region) level, respectively, for mG> = 1000 GeV. 
Remember, the iso-contours of the values for 1,2, and 3a assume the given model, while 
the (x^) chain value gives an indication of the overall fit. In this case, we get a somewhat 
worse fit to both experimental results than in the SM, with (x^)chain ~ 1-84. Fig. El^b) shows 
the estimated parameter space contours with an integrated luminosity of 10 fb~^, assuming 
the central values of both experimental measurements are not changed. The quality of 
the fit is marginally better than in the SM, with {x^)chain = 4.20 vs. 4.22 in the SM. We 
observe that the regions where fi or ql are small provide the best fit. The boundaries of 
all three contours can be understood from the theoretical predictions in Figs. M^c) and (d). 
To explain the discrepancy in the total cross section and in the asymmetry, values of fi and 
(?i in the top- left or bottom-right quadrants would be preferred. However, couplings here 
inevitably worsen the M^t distribution. In the top-right and bottom-left quadrants the fit 
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to the Mtt distribution is improved for small couplings , [(^lI ~ 1) but the agreement 
with the total cross section is slightly worse and an asymmetry smaller than the SM value 
is generated. For intermediate couplings in these two quadrants (1 ^ , Ig^l ^ 2), the 
fit to the Mft distribution is improved but the total cross section is reduced too much and 
the asymmetry is not improved significantly. Eventually, at large values of the couplings in 
these quadrants , IsfLl ^ 2.5), a large asymmetry is generated and the fit to the total 
cross section is improved but the distribution is greatly worsened. Furthermore, we note 
that the bands along the Ql = axis are slightly wider than those along the fL = axis due 
to the asymmetric contributions of and (^l to Tq', cf. Fig. |^e). 



2. mc = 2000 GeV 

When the G' boson is very heavy, only the interference term in Eq. f l2ip contributes to 
ApB, leading to 



A 



INT 
FB 



(/l + //?) {Ql + 9r)' 

This dependence is illustrated in Fig. [71 In order to get positive corrections to and 
the top pair production cross section, the product fiQi needs to be negative to compensate 



G' 



which is what we 



the negative sign from the denominator of the propagator — m; 
see in the upper left and lower right quadrants of Figs. [7] (a) and (b). The results of the 
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FIG. 7: Theoretical prediction contours for mc = 2000 GeV with /r = = in the plane of Jl 
and gL with fR = gR = 0: (a) cross section (pb) and (b) A^p*^. 
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FIG. 8: Same as Fig. [6] but for niQ' = 2000 GeV. The region bounded by the sohd curve defines 
the la region, the areas bounded by the dashed hues are within 2a, and the areas bounded by the 
dotted lines are within 3a. 

MCMC scan are shown in Fig. [8] for nic = 2000 GeV. The contours look quite different 
from those in Fig. O We note that values of and ql in the upper-left and lower-right 
quadrants are preferred which is where a large positive asymmetry is generated as seen in 
Fig. [71 This shows that the distribution is less constraining than in the m^/ = 1 TeV 
case as one would expect. The fit to the three experiments gives (x^)chain ~ ^-^^ 
current integrated luminosity which is slightly better than in the SM where (x^)chain ~ 1-75. 
The fit is improved relative to the SM at the upgraded luminosity with (x^)chain ~ 3.82 if 
the central values do not change as compared to the SM value of (x^)chain ~ 4.22. 

In Fig. [HI we observe that the upper-right and lower-left quadrants are not as tightly 
constrained as in the mc = 1 TeV case in Fig. [6] Here, the distribution is improved.^ 
In the case of the 1 TeV G' large couplings in these quadrants decrease the total cross section 
and asymmetry too much. However, in the 2 TeV case for large couplings in the upper-right 
and lower-left quadrants the NFS term becomes important due to the large width effects and 
this can mitigate the negative contribution to au and ApB from the INT term allowing for 
a better fit. This is why the 2 and 3a regions are not tightly constrained in the upper-right 
and lower-left quadrants for a 2 TeV G' in Fig. [8] Again, we note that the bands along the 

^ Naively, one might expect that the Mtf distribution would not be very important for a G' with mc = 
2 TeV. However, for couplings ^ 4 the width of the G' can be comparable to tog' (for = 

4, gL = 0, Tg' ^ O.Smc') and the 2 TeV G' can contribute to the 800 GeV < Mtt < 1400 GeV bin. 
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FIG. 9: Theoretical prediction contours for an axi-gluon of a^ps = ^NP (since ajNx = 0) and 
A% for (a,b) mc = 1000 GeV and (c,d) mc = 2000 GeV. 

Jl and ql axes are not symmetric due to their asymmetric contributions to Tg'- 



D. Axi-gluon: = -Jl and gn = -ql 

Now we study the axi-gluon case, in which G' only has axial couplings to the quark sector. 
This type of model has been explicitly proposed as an explanation of the ApB measurement 



lOj . There, the SM prediction for the 



without significantly affecting the total cross section 
ti cross section was taken to be larger than our value due to differences in and including 
incomplete NNLO calculations. Therefore, they did not need a significant correction to the 
cross section. 

In the axial limit, only the asymmetry-generating term of the INT in Eq. (!22|) remains. 
In general, all terms in the NPS remain. Therefore, at large mc'-, the INT term dominates 
and a rather large asymmetry can arise without a sizable contribution to cr^f or to At 
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lower values of mc/, the NPS term is increasingly relevant. 

1. mG' = 1000 GeV 

We show theoretical contours of cnps (since ajNT vanishes when integrated over cos^ 
in the axi-gluon case) and A°p^ in Fig. [9](a) and (b) for an axi-gluon of mass 1 TeV. We 
observe that a positive asymmetry is generated when the product fg, with f = Jl = —Jr 
and g = Ql = —Qr, is negative as we expect from Eq. [211 In Fig. [TUlfa) and (b), we 
perform MCMC scans and find a fit with (x^)cham = 1-56 for the current luminosity and 
{x^)chain = 2.94 for 10 fb~^ assuming the central values of the measurements do not change. 
Small values of either f or g are preferred due to the constraint on the Mu distribution. 
These values of (x^) chain for the 1 TeV axi-gluon are better than the corresponding values 
in the SM fit. This indicates that a 1 TeV axi-gluon has less tension with the current data 
than the SM and would also offer an improvement over the SM if the central values of the 
data remain the same with an upgraded luminosity of 10 fb~^. 

2. mc = 2000 GeV 

For a 2 TeV axi-gluon, we plot theoretical contours of cr^ps and ApB in Fig- M (c) and 
(d). In Figs. [TOT c) and (d), we show the results of a scan in the case of an axi-gluon with 
rriG' = 2000 GeV. The fit shows better agreement with the data in this case than in the SM, 
with {x^)chain = 1-15 for the current luminosity and {x^)chain = 2.04 for 10 fb~^ if the central 
values do not change. Due to the lessening of the M^f constraint for a heavier axi-gluon, the 
scans show somewhat different structure than in the 1 TeV case. The la allowed regions are 
located in the quadrants where fg<0 which is where a positive Aps is generated as seen 
in Fig. [9] (d). The regions of large coupling are constrained by the Mu distribution. There 
is again a slight asymmetry in the width of the allowed regions near the / and g axes due 
to the asymmetry in the width. 

Our results suggest that a heavy axi-gluon can offer a good explanation of the large App 
observed without increasing the disagreement in the M^f distribution too much, as proposed 
in Ref. M. 
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f f 

FIG. 10: Correlation of couplings for niQi = 1000 GeV (top panels) and mc = 2000 GeV (bottom 
panels) in the axi-gluon limit with the current integrated luminosity (left panels) and / Cdt = 
10fb~^ (right panels) assuming the central values of the measurements do not change. 

E. Other combinations of couplings 

Now let us study different combinations of couplings, e.g., /l = S'l = 0, /l = = and 
fR = gL = 0. Fig. [TT] shows the MCMC scan results of various combinations of couplings 
with rriG' = 2000 GeV and with the current luminosity. Purely right-handed couplings in 
the q-q-G' interaction give rise to the exactly same result as purely left-handed couplings, 
cf. Fig. IHa). But mixed combinations of left-handed and right-handed couplings, e.g., 
91 = /i? = and = /l = 0, result in a worse fit, with {x^) chain ~ 1-92, which is worse 
than the SM. This is mainly due to the INT effects which are sensitive to the signs of 
couplings; see Eqs. fl2T]) and f l27|) . Choosing /l = fi'R = or /r = Ql = causes the INT 
effects to generate a negative A^/^, leading to the bad fit. 
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IV. EFFECTIVE FIELD THEORY 



For a 2 TeV G' boson, due to the broad decay width, the NFS contribution is still sizable 
for large couplings in the above MCMC scans. It is interesting to ask what the effects are 
if only the INT term contributes. To that end, in this section we consider dim-6 effective 
operators that can interfere with the SM top quark pair production channel qq ^ g ^ ti. We 
further assume that the scale of the new physics is large enough that the NFS contributions 
(i.e. oc with A the new physics scale) can be neglected. For illustration we focus 

only on the operators which couple left-(right-)handed light quarks to left- (right-) handed 
top quarks so that contact with Sec. IIII Cl ean be made. They are listed as follows: 

Oti'^ = {qi,t\) {Q^n^Q) , (40) 

= (qi,t\^q) {Qrt\^Q) , (41) 

0if ^ = {u^.t^u) {trt^t) , (42) 

= [d^.t^d) {mH) , (43) 

where q and Q denote the SU{2)l doublets of the light (first two generation) quarks and 
heavy (third generation) quark, respectively, and u{d, t) are the right-handed gauge singlets. 
Here, and are the SU (3) and SU (2) matrices; appropriate contractions are understood. 
The first index in the superscripts of operators labels the color octet and the second index 
denotes the weak isospin. Other color and weak singlet operators are omitted as they cannot 
interfere with the SM channel. 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

/fi /fi /l 



FIG. 11: Correlation of couplings for mc/ = 2 TeV with the current integrated luminosity for (a) 
9L = fL = 0, (b) gji = fL = 0, (c) gL = fji = 0. 
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The effective Lagrangian of the four fermion interaction qqtt is thus given by 

q q 

^^'^^ = gl'^ (qi.PLq) ih^Pit) + g,'^ (qif^PRq) {h,PRt) , (44) 

where we exphcitly factor out a strong couphng strength g'^, and the reduced coefficients are 
given as follows: 

_ ^(8,1) , M8,3) d _ ^(8,1) _ ^(8,3) u _ MS) d _ MS) 

'^L — ^qq + ^qq ' '^L — ^ qq ^ qq ' '^R — ^ut ; I^R — ^ dt ■ 

Here the SU (3) generators are omitted and A denotes the new physics scale. 

The differential cross section of the EFT can be easily derived from Eq. ( l2Ti) by taking 
the limit of mc = A ^ 

f^L + '^R 



EFT _ 7r/3ag 

•^INT — j^g 



A2 



{(2-/32) +2/3 COS ^ + (/3 COS ^)^}. ^45^ 



Obviously, the coefficients n^iji only affect R but not A^^. We extract the cutoff scale and 
coefficients as follows: 



aeft _ + j Til^al 

(^)^ I 18(TeVf 



(2-/32) +2/3 cos e + (/3 cos ^) 2] I . ^4g^) 



which yields, after integration over s and convolution with PDFs, 

.eft _ ^ -i? ^^^^^ ^ ^^^^ ^ ^ ^^^^ ^^g2 ^) _ 



(tw) 



rfcos6' A \2 



The parameters {Aeft, Beft, Ceft) are listed in Table UTTl for various choices of factoriza- 
tion scale, where Ap^^ and aEFT are evaluated using Eq. fl26l) . It is clear that one needs 
positive kI or to get positive Afb but this inevitably gives rise to a negative contribution 
to a{tt). Hence, it is difficult to fit both the asymmetry and the top pair production cross 
section simultaneously at the la level. 

Since one cannot separate the coefficient k, from the cutoff A, we scan over the combination 
K/(A/TeV)2 and hmit ourselves to the region of \K/{A/TeVf\ < 10 in the MCMC scan. 
In Fig. [12] we plot the correlations between ^/(A/TeV)^ and 0"^* (top row) and between 
/t/(A/TeV)2 and A°pE (bottom row). For the current luminosity, the fit quality of the EFT 
is worse than the SM, (x^)chain ~ ^-80. The fit is marginally better than the SM, (x^)chain ~ 
4.20, for an integrated luminosity of 10 fb~^ if the central values of the measurements remain 
the same. The fit is worse than that of the 2 TeV left-handed G' due to the lack of a NFS term 
to balance the contributions of the INT term. This indicates the importance of resonance 
effects in the fit. 
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FIG. 12: (Left): Fitting contours in the EFT with the current integrated himinosity: (a) in the 
plane of A% and K/(A/TeV)2, (b) in the plane of a{tt) and ^/(A/TeV)^. (Right): same as the 
left column but for / Cdt = lOfb"^. 

V. FLAVOR-CONSERVING Z' BOSON 

An additional Z' can generate a nonzero Ai?^ if its coupling to quarks does not respect 
parity, 



Z'qq : le-f^ UlPl + IrPr) , 
Z'tt : iei^{gLPL + gRPR)- 



(47) 
(48) 



where e denotes the electromagnetic coupling strength. In contrast to G\ there is no in- 
terference between the Z'-mediated top quark pair production and the SM process. Even 
though the Z' amplitude interferes with the SM process qq — )■ 7*/Z* — tt, the latter contri- 
bution is negligible at the Tevatron. Only the NP resonance itself contributes to ApB when 
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TABLE III: Parameter {Aeft, Beft, Ceft) for EFT where A^/j^ and calculated after 

integrating over the angle 6. 



uu — )• tt 


dd tt 




Aeft 


Beft Ceft Ap^ 


a EFT 


Aeft 


Beft 


Ceft Af%^ 


a EFT 


mt/2 


-0.294 


-0.256 -0.092 0.395 


-0.648 


-0.052 


-0.040 


-0.040 0.355 


-0.113 


nit 


-0.215 


-0.185 -0.066 0.392 


-0.473 


-0.037 


-0.0288 


-0.10 0.353 


-0.062 


2mt 


-0.165 


-0.141 -0.050 0.389 


-0.363 


-0.028 


-0.022 


-0.007 0.350 


-0.082 



the collider energy is large enough to see the resonance effects. We consider the case where 
both the up and down quarks are gauged, but it is also possible to gauge the up and down 



quarks differently 



44|. 



Since the interference is absent for the color singlet Z\ only Anps contributes to NPS. 
For the s-channel diagram, the differential cross section of Z' can be easily derived from 
that of G' by omitting the color factor 2/9 in Eq. fl22|) and replacing by aemi yielding 



da 



dcos9 



7r^9 
7r/3a 



2 

em 



fi2 



If + mm 



9l 



C- G 



X 



9l 



9l 



(1-/3^) 



Xfl-fl) (gl-gl) 



(3cose + {(3cosef\ . (49) 



'(/K/I) i9l + 9iy 

Negative searches for the Z' boson at the Tevatron impose several lower bounds on the Z' 
mass, roughly above 1 TeV for couplings of order electroweak size. For a leptophobic Z' 
boson, the bound is slightly looser, mz' > 700 GeV. Owing to the rapid drop of the PDFs, 
the Z' boson contributes significantly only in the resonance region, where /3 — > 1. Further 
noting that the coefficient of the (1 — term in Eq. is always less than one, we can 
drop this term and obtain 



oc 



ifl - Pr) (gl - gl) 



^ ifl + Pr) (gl + gl) ' 

As in the G' study, we turn off the right-handed couplings first. Fig. 



(50) 



displays the 

correlation between and ql couplings for a 1000 GeV Z' boson with fR = gR = 0. Like 
the 1 TeV G', the distribution favors smaller or ql. The fit is worse than in the 



SM with (x^) 



chain 



1.85 for the current luminosity and (x^) 



chain 



4.28 for 10 fb"^ if the 



central values remain the same. For a 2 TeV Z' boson with = gR = 0, we see the results 
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FIG. 14: Same as Fig. [6] but for 2 TeV Z' . 

of the MCMC scan in Fig. [H The fit is somewhat better than the SM: ^od) chain = 1-62 
for the current luminosity and {'}d)chain ~ '^■^^ f*-*^ lessening of the 

importance of Mti for the higher mass Z' . Furthermore, large couplings are allowed at the 
1(T level. However, note that the unitarity constraint derived for the process uu ^ Z' ti 
requires ^ 28; see Appendix |A] for further details. Now, the heavy Z' contributions are 
very sensitive to width effects. For Jl = Ol = 10 we obtain Tz' — O.SSm^/. The positive 
contributions to the Mu distribution, particularly in its last bin, are somewhat constrained. 
We conclude that this model can offer a small improvement over the SM in describing A^^, 
ajf, and the tt invariant mass distribution simultaneously. 
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VI. FLAVOR- VIOLATING Z' AND W'^ MODELS 



In this section we consider a flavor- violating Z' model, which includes a u-t-Z' interaction. 
Such a FCNC could appear at tree-level or loop-level. Rather than focus on a speciflc model, 
we consider the following effective coupling of u-t-Z' jsj: 



C = cu^UlPl + fRPR)tZ'^, 



(51) 



where e denotes the electromagnetic coupling strength. In addition to the SM QCD pro- 
duction channel, uu ^ g ti, the top quark pair can also be produced via the process 
uu — 7- ti with a t-channel Z' boson propagator. The top quark asymmetry is naturally 
generated by this new process which also interferes with the SM production mode. There- 
fore, the differential cross section versus the cosine of the top production angle is 9 given as 



follows 



9|: 



da 



dcos9 



•Asm + AfrJT + A^ 



where 



INT 



A'' 



NP 



727rs s{t-Ml) 
1287rs {t - Ml,y 



2(u 



m 



2\2 



NPi 



m 



+ 2smf + -—^{{t - mfy + sm 



Ml 



(52) 



(53) 



{al + gV){u - miY + 2gigj,s{s - 2m 



{gl + glf{UMl, + {t-Ml,] 



(54) 



and Asm is given in Eq. f l20|) . The interference between the QCD and EW processes can be 
easily understood as follows. The SU (N) gluon propagator can be split into a U{N) gluon 
propagator and a U{1) gluon propagator 45l |. 



(55) 



where the U{1) gluon, carrying a factor is unphysical. Color flow of the SM QCD 

channel (i.e., u ^ t and t — n) is then exactly the same as the Z' induced t-channel 
diagram, resulting in interference between both processes. 

Within the SM, the FCNC coupling u-t-Z vanishes at tree-level, but can be generated 
at one loop. However, the one-loop generated coupling is strongly suppressed by the GIM 
mechanism, making the FCNC top interactions very small. In models beyond the SM this 
GIM suppression can be relaxed, and one-loop diagrams mediated by new particles may 
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also contribute, yielding effective couplings orders of magnitude larger than those of the 
SM. Since the coupling strength of this FCNC interaction is typically at the order of the 
SM weak interaction, the coefficients fi and fji are expected to be much smaller than 1. 
Therefore, it is not easy to generate a large asymmetry from a loop-induced u-t-Z interaction. 
However, the couplings and could be larger if they are generated at tree-level. 

While the value of fn is not well constrained by direct or indirect search experiments, 
the value of fi is tightly bounded by the _B-sector. The left-handed coupling in Eq. f lSTj) 
originates from the gauge interaction of the Z' boson to the first and third generation quark 
doublets, 

£ = Qlii^D^Ql + h.c. , (56) 

where qhiQi) denotes the first (third) generation quark doublet and the covariant derivative 
is = id^ + iB'Z'^, where B' is the charge. The flavor violating interaction di-hi-Z' then 
follows directly from the gauge invariance, which can contribute to the B^-B^ mixing at the 
tree level ^ . A coupling of the form 



C = edY{fLV:,VttPL)bZ'^ , (57) 

follows from Eq. [56] after rotating to the mass eigenstate basis with fi as in Eq. [51] and Vud 
and Vtb elements of the CKM matrix. Assuming no additional NP effects arise, this gives a 
contribution to the mass difference between B^-B^ of 

where fs^ is the Bd decay constant and B is the "bag parameter" that characterizes the 
deviation from the vacuum saturation approximation. If we conservatively require that this 



contribution does not exceed the experimental value of 3.34 x 10 MeV 47|], we can set a 
limit on of 



where we use fsivB = 216 ± 15 MeV [48|] and I^^J^Vtb] ~ 1 [47|. As a result, we choose 
fL = hereafter. Furthermore, the unitarity constraint derived for the process ut —> Z' ^ ut 
requires only ^ 28; see Appendix Rl for further details. 



2 



A similar correlation among the gauge boson and the third generation quarks in the SM has been studied 
in Ref. [461. 
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FIG. 15: Value of the same-sign top quark production cross section at the Tevatron for a flavor 
violating Z' with = 0, //j = 1 (Eq. (f5T]l ). 



The most striking signal of the FCNC Z' model is same-sign top pair production through 
the processes uu — tt via a t-channel diagram mediated by the Z' boson. Recently, the CDF 
collaboration searched for the same-sign top pair signature induced by the maximally flavor- 



violating scalars at the Tevatron and found no evidence of new physics beyond the SM [49 1. 
In their analysis, the upper limit to the production cross section of same-sign top pairs is of 
the order of 0.7 pb. We show, in Fig. [151 the same-sign top production cross section at the 
Tevatron for couplings //j = 1 and /l = 0. The cross section scales with the right coupling 
as cr(tt + tt) ~ if //^ = 0. Direct production via t-channel Z' exchange dominates and 
is severely constrained as the couplings increase. Note that we do not consider the possible 
effects that same-sign top production could have on a measurement of Aps which requires 
a delicate analysis and will be presented elsewhere [50|]. 

In Fig. [16] we show the result of the MCMC scan in the plane of mz' and (a) for 
the current integrated luminosity and (b) for expected 10 fb~^. For the current integrated 
luminosity, we impose the constraint of a{tt + it) < 0.7 pb 49|, whereas for 10 fb~^, we 
assume the cross section limit scales with l/y/C, giving a{tt + tt) < 0.4 pb. The value of 
(x^) chain ~ ^'^^ curreut integrated luminosity indicates that the FCNC Z' model fits 

a{ti) ApB, and the Mu distribution worse than the SM. Note the quality of fit is maintained 
even if we fix the Z' mass to be specific values as in the flavor-conserving G' and Z' cases. 
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FIG. 16: Correlation of fn and mz' with = 0. The left panel is for the current integrated 
luminosity while the right panel is for an upgraded luminosity of / Cdt = lOfb"^. Tension with 
the same-sign top production constraint from the Tevatron prevents this model from becoming a 
better fit than the SM. 

Since the INT effects lead to a negative asymmetry, one needs a large NPS contribution to 
overcome the negative INT contributions to generate the positive asymmetry. That requires 
a very large coupling, as seen in the upper and lower la contours of Fig. IT^ a) and (b), 
which is near the constraint of a{tt-\-ii) < 0.7 pb. In this model, the predicted value for the 
same-sign top pair production cross section is pushed to just below the limit taken. Overall, 
while there is tension between the positive asymmetry and small a(tt + it), we find a fit 
not much worse than the SM. With 10 fb~^, the fit remains worse than in the SM with 



chain 



4.39. 



The observed top asymmetry may also be induced by a flavor-changing interaction via 
a charged W boson |9|, [l7| . The top quark pair can be produced in the channel dd — )■ ti 
via a t-channel W boson propagator. As in the flavor-violating Z' case, we consider the 
following effective d-t-W coupling: 



C = edY{fLPL + fRPR)tW' 



(60) 



where e denotes the electromagnetic coupling strength. The differential cross section of 
dd ^ ti is the same as Eq. (15^ with the substitution u{u) -> d{d). 

One advantage of the flavor-violating W model is that it does not suffer from the con- 
straint of same-sign top pair production at the Tevatron. Fig. [17] displays the correlation 
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FIG. 17: Correlation of couplings for m^y = 2000 GeV. The left panel is for the current integrated 
luminosity while the right panel is for an upgraded luminosity of J Cdt = lOfb"^. The la region 
is bounded by the solid curves, the 2a by the dashed curves, and the 3a by the dotted curves. 
Couplings of either /l ~ and /r ~ ±10 or and Jl ~ ±10 are consistent with the Aps 

and att measurements. 

between Jl and fn couplings for a 2000 GeV W boson. For the current luminosity the la 
contours are symmetric for and fn and the innermost one includes fL = fR = 0- In order 
to generate positive asymmetry, the couplings and fji need to be large enough to over- 
come the negative INT contributions. The typical values of couplings fi and //? are in the 
range of ±10 to ±20. For the current luminosity, {x'^)chain = 1-67 which indicates a better 
overall fit than the FCNC Z' boson due to the lack of the same-sign top constraint, and a 
better fit than the SM. Due to the PDF dependence, the (i-quark initiated tt production 
via the W is smaller than the w-quark initiated production through a Z'. Therefore, larger 
couplings are required to maintain the production cross section than in the Z' case. With 
upgraded luminosity, {x^) chain ~ "^-^^ provided that the central values of ^"^^ '^tf 
maintained which offers more improvement over the SM. 

We focus on a heavy W due to general constraints from electroweak precision and flavor 
measurements. In general, a W is associated with a broken non-abelian gauge group and 
one must also consider a neutral gauge boson, Z', whose mass is typically degenerate or 
nearly so with that of the W. If this Z' has predominantly flavor- changing couplings to 
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top quarks, then it falls into the previous case we analyzed. If its coupling to top quarks 
is flavor-conserving then one would expect to produce top quark pairs through s-channel 
Z' exchange. However, such a process suffers from PDF suppression and is negligible in 
comparison to the W contribution considered above and therefore we ignore it here. ^ 



VII. FLAVOR- VIOLATING SCALAR SjS^ 

In addition to spin-1 exchange, we also consider the FCNC top interaction with a new 
color singlet scalar S': 

L D eS' ihuntL + jRULtR) , (61) 

where S' is an SU (2) doublet and we parameterize the overall coupling strength with respect 
to the weak coupling e. If we assume S' to be the SM Higgs boson, then the FCNC top 
interaction originates from the dimension-6 operator 

_ 

Results for this operator can be obtained from those for Eq. flMl) with the substitution 

2 

e — > As will be shown below, such an effective coupling is too small to generate a sizable 
asymmetry however. Hence, it is difficult to explain the asymmetry with the SM Higgs 
boson effective coupling without introducing additional heavy scalars. 
The differential cross section is written as 



da 



dcos9 



A. 



SM 



A 



S' 

NPS 



(62) 



where 
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A^' 
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is t 



27ra,aem/3 (/| + fl) 


smf + 
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vr«L/3 ifl + fir 




-m|,) 
[t - mf 


f 



(63) 
(64) 



^ A model with a light W' and Z' has been proposed in Ref. [17| and may lead to a naturally good fit. 
This model has potential implications for precision electroweak observables which have not yet been fully 
explored. 
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with = e^/(47r). Due to the repulsive scalar interaction, the NPS contributions generate 
a negative Ap^. In order to generate a positive Ap^, the scalar S' needs to be very light, 
generally ms' < mt, and to have large couplings with the top quark. However, such a light 
scalar leads to a new top quark decay channel t — S'u, which is tightly constrained 47 |. 
Therefore, we consider scalar masses that are larger than top quark mass to forbid this 
new decay channel. Furthermore, the flavor-violating coupling will be highly constrained 
by D^-D^ mixing (a AC = 2 process) if one assumes a universal flavor-violating coupling 
among the three families of quarks. However, from a purely phenomenological perspective, 
we assume that the second generation quarks are not involved in the flavor- violating Yukawa 
interaction which leads to no constraint on the Yukawa couplings Jl and In other words, 
fut (/L;/ij) is taken as a free parameter and is only constrained by considerations of unitarity 
(see Appendix Rl for details). 



We plot Ap^ as a function oi f = \J f1 + for a scalar of mass 2 TeV in Fig. [18] (d) and 
see that it is indeed negative. This does not pose a problem with respect to the measurement 
of a positive asymmetry since the scalar interferes destructively with the SM which gives a 
negative i?, defined in Eq. (I6l) and shown in Fig. [18] (c). Thus, the total asymmetry, which 
is related to A^^ and R in Eq. (|5|), is positive. This is plotted in Fig. [18] (a). 

We also consider a charged scalar, S*"*^, which couples to top quarks as in Eq. (!6T|) . but 
with the replacement u ^ d. In Fig. [TQ] we see the result of the MCMC scan in the plane of 
/l and /r for a 1000 GeV neutral 5" and charged S"^ scalar: (a,c) for the current luminosity 
and (b,d) for expected 10 fb~^, again assuming the central value of experimental data is not 
changed. For the neutral scalar, we impose the same-sign top pair production constraint 
discussed in Sect. [VT] We plot the same-sign top pair production cross section for a neutral 
scalar with = //? = 1 at the Tevatron in Fig. [20] In the case of the neutral scalar, with 
the current luminosity we find a marg inally better fit than in the SM, (x^)cham = 
The fit improves relative to the SM with (jk^) chain ~ 3-15 for 10 fb~^ if the central values 
remain the same. For the charged scalar we obtain a similar fit at the current luminosity 
with (x^)c/iam ~ 1-72. At 10 fb~^ with the same central values we obtain a slightly better 
fit) {'X^) chain ~ 3.04, than in the case of the neutral scalar due to differences between the u 
and d PDFs. Although the negative (J/att decreases the total cross section ajf, it allows for 
good agreement with the distribution and gives a positive contribution to which 
results in a better overall fit. 
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FIG. 18: Theoretical expressions for a flavor-violating neutral scalar S. (a) A^p^^ (cf. Eqs. ©-([S])) as 
a function of / = \J~ff^^r~fji fo^' a flavor violating scalar with ms = 2 TeV. (b) a^p = ctjnt + ctnps 
vs. /. Note that the interference is destructive, (c) R, as deflned in Eq. ([6]). Note that it is negative 
for all /. (d) Apj^ vs. / which is negative, as mentioned in the text, (e) The ratio of (Tnps to 
(Tint- The INT term dominates over the NPS term for a scalar of mass 2 TeV. 

VIII. CONCLUSION 

We have examined a number of models for new physics in top quark pair production 
which could account for the larger-than-expected forward-backward asymmetry observed at 
the Fermilab Tevatron, while not significantly disturbing the approximate agreement of the 
cross section au and its distribution with the Standard Model predictions. Our results 
are summarized in Table IIVI 

The results summarized in Table IIVI show that it is not easy to account for the larger- 
than-expected value of ApBlti) observed at the Fermilab Tevatron while maintaining the 
good agreement between theory and experiment for the production cross section au and 
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FIG. 19: Correlation of couplings for (a,b) ms' = 2000 GeV and (c,d) ms± = 2000 GeV. The 
left panels are for the current integrated luminosity while the right panels are for an upgraded 
luminosity of / Cdt = 10 fb^^. The regions of each figure, from innermost to outermost, are (a) 
within la, within 2a, within 3a, and greater than 3a; (b) within 2a, within la, within 2a, within 
3a, and greater than 3a; (c) within la, within 2a, within 3a, and greater than 3a; (d) within 2a, 
within la, within 2a, within 3a, and greater than 3a. The couplings are varied only within their 
allowed values from unitarity considerations (see Appendix [A|) . 

differential rate 

dMtt 

Of the models considered, those that provide a fit better than the SM for the applicable 
data are a 1 TeV or 2 TeV flavor-conserving G' with axial couplings, a 2 TeV W (or a 2 TeV 
flavor- conserving G' or Z') with chiral couplings. Other models we considered provide at 
most a mild improvement with respect to the SM case. The 1 TeV cases often give large 
corrections to the M^^ distribution since the additional signal is well inside the data region. 
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FIG. 20: Value of the same-sign top quark production cross section at the Tevatron for a flavor 
violating S' with = fB = ^ (Eq. (j6T]) ). The gray region indicates ms' < mt. 

Finally, in Table |Vl we examine in detail the contribution to the {'X^)chain in the small 
neighborhood around the best points in parameter space of the axial G', chiral G', FV chiral 
W\ and FC chiral Z' models as well as the SM. Assuming the current integrated luminosity 
for each measurement and a mass of the new states responsible for the ApB deviation of 2 
TeV, we find: 

• The 2 TeV axial G' model provides the best overall fit to the experimental data 
considered in this article. It improves the agreement with experiment of both the 
total top quark production cross section and the forward-backward asymmetry. The 
fit to the Mtt distribution is slightly worse than in the SM case but still in good 
agreement with data. 

• The 2 TeV W can also lead to a good fit to the data. It is able to generate a large 
asymmetry and to improve the agreement of total cross section with data without 
disturbing the differential cross section sizably for some regions of parameter space. 
However, large couplings are needed. Note that in Table |Vl the best-fit point in the 
W case has a lower than the axial G' although the {y^) chain is lower for the axial 
G' indicating that the axial G' gives a better overall fit. Stated differently, the W 
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TABLE IV: Models for top quark pair production and their ability to account simultaneously for 
the cross section and forward-backward asymmetry in top quark pair production at the Tevatron. 
Flavor Conserving (FC) and Flavor Violating (FV) models are considered. 



Model Result 

FC G' 1 TeV Poor fit for fR = gR = considered due to constraint; 

Good fit for axial couplings 
2 TeV Fit not significantly improved with respect to the SM for fR = gR = 0; 

Excellent fit for axial couplings 
EFT Poor fit; ApB consistently smaller than measured value. 

FC Z' 1 TeV Poor fit due to constraint 

2 TeV Fit not significantly improved with respect to the SM for fR = gR = 
FV Z' Poor fit due to ^jj— constraint and same-sign top constraint 

FV W 2 TeV Good fit although large couplings necessary with a large amount of fine-tuning 
FV S", 5^ 2 TeV Tension with small predicted a^j leads to poor fit with current data; 

however, a good fit would be obtained if central values were unchanged after 10 fb 



requires a greater amount of fine-tuning of its parameters to fit the data than the 
axial G'. This is seen in the large value of 5 {y^) point'-, a slight perturbation of the 
best fit points greatly decreases the quality of the fit. The W model provides such 
a large value since the contributions from (Jti and Aps are aligned and increase 
together with couplings that deviate from the minimum couplings. This is to be 
contrasted with the other models in which the increasing contribution from, say, 
ApB ^^ compensated by a smaller contribution from ati-, resulting in a total value 
that remains relatively flat. 

• The 2 TeV chiral G' and Z' do not lead to significant improvement over the SM. They 
reduce the discrepancy with the asymmetry measurement although they are unable to 
reduce it below 2(7 without disturbing the Mti distribution due to their large widths. 

• The 2 TeV FV scalars 5" and have fits that are not significantly improved with 
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TABLE V: Values of for selected models. Assuming the current integrated luminosities, the 
{'X^)point contributions from the production cross section, forward-backward asymmetry and in- 
variant mass distribution in top quark pair production at the Tevatron are given for the chain only 
in the small neighborhood of the points that provide the best fit to the data in each model. The 
value of {x^)point is the Total/ Ndof, where N^of = 3, and provides a measure of how well the best 
fit point in each model fits the data. The value 5{x'^)point is the {x'^)point value of a box that is 
±10% wide in the couplings for the best fit point and is a rough measure of fine tuning. Note the 
SM value for (x^)pomt is equivalent to that given in Eq. [16l The (x^)chmn values for each model 
are also listed. 



Model 




Afb 


da 
dMtt 


Total (x^)pomt 


^{X^) point 


{X } chain 


FV W ih -- 


= 17.0, /i? = 0.0, miy' = 


= 2 TeV): 


0.63 


0.06 


0.27 


0.96 


0.32 


41.4% 


1.67 


axial G' (/ 


= 2.5,5 = -2.5,mG' = 


-- 2 TeV): 


0.75 


0.47 


0.75 


1.97 


0.66 


1.6% 


1.15 


chiral G' (/l 


= 2.0,(7L = -2.0, 


= 2 TeV): 


0.59 


2.91 


0.62 


4.12 


1.37 


0.4% 


1.69 


EC Z' ih 


= 8.0,(7L = 8.0,mz' = 


2 TeV): 


0.67 


2.66 


0.91 


4.24 


1.41 


1.0% 


1.62 


FV S' Ul 


= 7.0, fR = 0.0, ms> = 


2 TeV): 


2.01 


2.85 


0.09 


4.96 


1.65 


0.2% 


1.72 


FV 5± {fL 


= 7.5jR = 0.0,ms' = 


2 TeV): 


2.01 


2.85 


0.12 


4.98 


1.66 


0.2% 


1.72 


SM : 


1.12 


4.07 


0.06 


5.25 


1.75 




1.75 



respect to the SM. They lead to a significant discrepancy in au and only slightly 
improve the fit to ApB and da/dM^ with respect to the SM. 

In this work, we have used the full NLO QCD ti production cross section. It is worth 



noting t 
in Ref. 



lat partial NNLO QCD corrections to the tt cross section have been calculated 
23| and give rise to an enhancement of the total cross section of about 0.3 pb. 
This indicates that higher order QCD corrections might improve the agreement between the 
measured total cross section and its value in the SM, and therefore areas of NP parameter 
space which give negative contributions to the total cross section will be less constrained. 
Such negative contributions, however, may be in tension with the observed invariant 
mass distribution. A detailed collider simulation, including the complete NNLO corrections, 
would be therefore highly desirable in order to make a more reliable comparison of the 
predictions of these models with data. 

41 



Crucial to the test of any model is also the accumulation of more integrated luminosity 
at the Tevatron, in order to demonstrate deviations from the Standard Model exceeding 3cr. 
Until then, the observed Aps values cannot be regarded as anything more than a hint of 
new physics. 
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Appendix A: Unitarity constraints 

In this appendix, we explore the unitarity constraints on new physics models considered 
in this work. The weak isospin amplitude Ai^ {I being isospin index) can be decomposed 
with respect to orbital angular momentum according to 

oo 

= 167r ^(2£ + l)PKcos 9)al (Al) 



With the normalization 53a| = |a{P, the unitarity constraint requires 

\^ai\ < ^, (A2) 

where a| could be projected via: 

4 = — r dcos9Pe{cos9)M^. (A3) 
327r J~i 

First consider the flavor-conserving Z' models, which involve the s-channel diagram only. 
Note that the constraints of the G' model can be easily derived from flavor-conserving Z' 
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model. The helicity amplitudes for qq — > tt are represented by A{Xq, \q, Xt, Xt\), where 
At = — ,+, respectively, indicates a left-handed and a right-handed top quark. Apart from 
the common factor 

s — vn^i 

the non- vanishing helicity amplitudes from the diagram mediated by the Z' boson are 





+, 




= /l sin e mt [ql + 






(A4) 


^(-, 


+, 




= -/l(1 + cos^)E[(1 + 




- (1 - l^t)gR] , 


(A5) 


^(-, 


+, 


+,-) 


= ^(l-cos^)E[(l- 






(A6) 
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+.+) 


= -/l sin e mt [ql + Qr^ 






(A7) 








= fn sin 6* rrit [gi + 5-/?] 






(A8) 








= fR{l-cose)E[{l + 




- (1 - , 


(A9) 








= -/ij(l + COS^)E[(l- 




h(l + A)^7/?], 


(AlO) 






+.+) 


= -/ij sin e mt [ql + S-R 






(All) 



where fit = y 1 ~ m^/ E"^. In the cm. frame of the tt pair the 4-momenta of the particles 
are chosen to be 

Pq = E (1, 0, 0, 1) (A12) 
Pq = E (1, 0, 0, - 1) (A13) 
pt = E{1, Asin^,0, Acos^) (A14) 
Pt = E (1,-Asin^,0,-Acos^). (A15) 

In the high energy limit E ^ mt, only (1 + Pt) terms contribute. The J = 1 partial-wave of 
A{—, — , +) amplitude is 

ai=\qq ^ Z' ^ tt) = ^hgL, (A16) 

yielding the following limit U'lOlI < Q/<^em- Similarly, one can derive the following con- 
straints 

|/l|<28 and < 28 

from the qq ^ Z' ^ qq and tt ^ Z' ^ ti processes. 
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Now consider the flavor-violating Z' model. We consider the scattering ui Z' ^ tu, 
the calculation of which is identical to the flavor-conserving Z' model. In the high energy 
limit ^/s ^ mt, we obtain the following unitarity bound on from the helicity amplitude 
A{+, -, +, -) in Eq.lAlOl 

Ifnl ^ 28 . (A17) 

Finally, we consider the scattering ut S ut to derive the unitarity constraint for the 
flavor- violating S model. The helicity amplitudes are represented by A{Xg, Af, Xu, At). In the 
high energy limit s ^ rrit, the non- vanishing helicity amplitudes are 

+,+,+) = e^fl (A18) 

A{+, +,-,-) = e^fnh, (A19) 

A{-, -,+,+) = e'hfR, (A20) 

^(-,- -,-) = e'fl (A21) 

where fi and are given in Eq. [61] The J = partial-wave of A{+, +, +, +) is 

ai=\ui^S^ut) = ^fn, (A22) 
yielding the unitarity limit /| < 2/aem- Hence, ^ 16, \fR\ ^ 16, and < 256. 



[1] CDF, G. Strycker et at, (2009), Public-Note:: CDF/ANAL/TOP/PUBLIC/9724. 

[2] CDF, T. Aaltonen et al, Phys. Rev. Lett. 101, 202001 (2008), arXiv:0806.2472. 

[3] O. Antunano, J. H. Kiihn, and G. Rodrigo, Phys. Rev. D77, 014003 (2008), arXiv:0709.1652. 

[4] DO, V. M. Abazov et al, Phys. Rev. Lett. 100, 142002 (2008), arXiv:0712.0851. 

[5] J. H. Kiihn and G. Rodrigo, Phys. Rev. Lett. 81, 49 (1998), arXiv:hep-ph/9802268. 

[6] J. H. Kiihn and G. Rodrigo, Phys. Rev. D59, 054017 (1999), arXiv:hep-ph/9807420. 

[7] A. Djouadi, G. Moreau, F. Richard, and R. K. Singh, (2009), arXiv:0906.0604. 

[8] S. Jung, H. Murayama, A. Pierce, and J. D. Wehs, Phys. Rev. D81, 015004 (2010), 
arXiv:0907.4112. 

[9] K. Cheung, W.-Y. Keung, and T.-C. Yuan, Phys. Lett. B682, 287 (2009), arXiv:0908.2589. 

[10] P. H. Frampton, J. Shu, and K. Wang, Phys. Lett. B683, 294 (2010), arXiv:0911.2955. 

[11] J. Shu, T. M. P. Tait, and K. Wang, (2009), arXiv:0911.3237. 

44 



[12] A. Arhrib, R. Benbrik, and C.-H. Chen, (2009), arXiv:0911.4875. 
[13] P. Ferrario and G. Rodrigo, (2009), arXiv:0912.0687. 

[14] I. Dorsner, S. Fajfer, J. F. Kamenik, and N. Kosnik, (2009), arXiv:0912.0972. 
[15] D.-W. Jung, P. Ko, J. S. Lee, and S.-h. Nam, (2009), arXiv:0912.1105. 
[16] J. Cao, Z. Heng, L. Wu, and J. M. Yang, (2009), arXiv:0912.1447. 
[17] V. Barger, W.-Y. Keung, and C.-T. Yu, (2010), arXiv:1002.1048. 



[18] CDF public note, http: //www-cdf . fnal.gov/physics/iiew/top/2009/xsection/TopXsNN/topXs_46fb_pul 



[19] CDF pubhc note, http : //www- cdf . f nal . gov/physics/new/top/ conf Notes/cdf 9913_ttbarxs4invfb .ps 



[20] P. M. Nadolsky et ai, Phys. Rev. D78, 013004 (2008), arXiv:0802.0007. 
[21] U. Langenfeld, S. Moch, and P. Uwer, Phys. Rev. D80, 054009 (2009), arXiv:0906.5273. 
[22] V. Barger, W.-Y. Keung, and G. Shaughnessy, Phys. Rev. D78, 056007 (2008), 
arXiv:0806.1962. 

[23] N. Kidonakis and R. Vogt, Phys. Rev. D78, 074005 (2008), arXiv:0805.3844. 
[24] P. Nason, S. Dawson, and R. K. Elhs, Nucl. Phys. B303, 607 (1988). 
[25] W. Beenakker, H. Kuijf, W. L. van Neerven, and J. Smith, Phys. Rev. D40, 54 (1989). 
[26] M. Cacciari, S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, JHEP 09, 127 (2008), 
arXiv:0804.2800. 

[27] S. Moch and P. Uwer, Nucl. Phys. Proc. Suppl. 183, 75 (2008), arXiv:0807.2794. 
[28] CDF, T. Aaltonen et al, Phys. Rev. Lett. 102, 222003 (2009), arXiv:0903.2850. 



[29] DO pubhc note, http : //www-dO . f nal . gov/RTan2Physics/WWW/results/prelini/TaP/T83/T83 . pdf 



[30] J. Pumplin et al, JHEP 07, 012 (2002), arXiv:hep-ph/0201195. 

[31] J. C. Rati and A. Salam, Phys. Lett. B58, 333 (1975). 

[32] L. J. HaU and A. E. Nelson, Phys. Lett. B153, 430 (1985). 

[33] P. H. Frampton and S. L. Glashow, Phys. Lett. B190, 157 (1987). 

[34] P. H. Frampton and S. L. Glashow, Phys. Rev. Lett. 58, 2168 (1987). 

[35] J. Bagger, C. Schmidt, and S. King, Phys. Rev. D37, 1188 (1988). 

[36] F. Cuypers, Z. Phys. C48, 639 (1990). 

[37] Q.-H. Cao, S.-L. Chen, E. Ma, and G. Rajasekaran, Phys. Rev. D73, 015009 (2006), 

45 



arXiv:hep-ph/0511151. 

[38] C. D. Carone, J. Erlich, and M. Sher, Phys. Rev. D78, 015001 (2008), arXiv:0802.3702. 
[39] P. Ferrario and G. Rodrigo, Phys. Rev. D78, 094018 (2008), arXiv:0809.3354. 
[40] M. V. Martynov and A. D. Smirnov, (2009), arXiv:0906.4525. 

[41] R. Contino, T. Kramer, M. Son, and R. Sundrum, JHEP 05, 074 (2007), arXiv:hep- 
ph/0612180. 

[42] J. M. Campbell and R. K. Ellis, Phys. Rev. D62, 114012 (2000), arXiv:hep-ph/0006304. 
[43] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, (2009), arXiv:0912.3375. 
[44] J. L. Rosner, Phys. Lett. B387, 113 (1996), arXiv:hep-ph/9607207. 

[45] F. Maltoni, K. Paul, T. Stelzer, and S. Willenbrock, Phys. Rev. D67, 014026 (2003), 

arXiv:hep-ph/0209271. 
[46] E. L. Berger, Q.-H. Cao, and I. Low, (2009), arXiv:0907.2191. 
[47] Particle Data Group, C. Amsler et al, Phys. Lett. B667, 1 (2008). 

[48] E. Gamiz, C. T. H. Davies, G. P. Lepage, J. Shigemitsu, and M. Wingate, Phys. Rev. D80, 

014503 (2009), arXiv:0902.1815. 
[49] CDF, T. Aaltonen et al, Phys. Rev. Lett. 102, 041801 (2009), arXiv:0809.4903. 
[50] Q.-H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy, and C. E. M. Wagner, in preparation. 



46 



